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ABSTRACT 


A time- dependent: computer model of radiative-convective- 
conductive heat transfer in the Martian ground- atmosphere 
(G-A) system was refined by incorporating an intermediate line 
strength CO 2 band absorption which together with the strong- 
and weak-line approximation closely simulated the radiative 
transmission through a vertically inhomogeneous stratification. 
About 33,000 CO 2 lines were processed to cover the spectral 
range of solar and planetary radiation. Effective absorption 
by silicate dust particulates, with height-dependent poly- 
disperse concentration functions , was built into the model to 
study the impact of dust on the ground-atmosphere temperature 
field as a function of time. 

This model was subsequently attuned to IRIS, IR-radio- 
metric and S-band occultation data for representative dates 
during the time period November 1971 through February 1972. 

Very satisfactory simulations of the measured Iris spectra 
were accomplished for the dust- free condition. In the case 
of variable dust-loads, the simulations turned out to be 
sufficiently fair such that some inferences into the effect 
of dust on temperature were justified. 

The numerical model simulated at high spatial, and 
temporal resolution (52 atmospheric and 30 subsurface levels 5 
with time-steps between 4 and 10 min) the heat transports in 
the G-A system. The algorithm was based on the solution of 
the radiative transfer and enthalpy rate equations. Ground 
and atmosphere were coupled by an internal thermal boundary 
condition. 

Basic findings were that high-resolution temperature 
soundings and their diurnal histories, generated by the 
attuned model, serve to adjust IRIS-in verted and S-band 
occultation-inverted temperature soundings together with IR- 
radiometric data. 
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1.0 INTRODUCTION 


In a previous study , the feasibility of reinterpreting 
Mariner 9 IRIS data by a radiative-convective- conductive 
transfer (R.C.C.T.) model was demonstrated for the case of 
a dust-laden Martian atmosphere (Pallmann, 1974). The 
findings of this modeling effort stimulated the proposal of 
an extended investigation into attuning the R.C.C.T. model 
to several modes of Mariner 9 data such as the IR-radiometric 
and the S-band occultation data together with the IRIS 
measurements . 

Various refinements of the CO ^ -transmission function 
were envisioned, attunable to the prevailing atmospheric 
conditions on Mars. Based on recent developments of the 
computational treatment of radiative transfer in Mie-seat- 
tering and absorbing planetary atmospheres (Canosa and 
Penafiel, 1973: Dave and Canosa, 1974: Dave, 1975) an Im- 
proved algorithm was incorporated into the R.C.C.T. model in 
order to represent the effective absorption of radiative 
energy by the silicate particulates during the major dust 
storm of 1971 (Martin, 1974). Various components of the 
dust algorithm were quantitatively adjusted to the findings 
of Conrath (1974) and the simulations of polydispersions 
offered by Deirmendjian (1969). 

An additional computer program, labeled "FLUXQUT" , 
generated the numerical values of the upwardly directed 
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hemispherical flux, at the top of the atmosphere, of the 
thermal emissions from the ground surface and all atmos- 
pheric layers, after due transmission through the inter- 
vening layers, within the spectral range 200-2,000 cm 

A particular subroutine for convective heat trans- 
port was incorporated into the R.C.C.T. model with the pro- 
vision that wherever the vertical transfer of (solar and 
planetary) radiative energy generated a superadiabatic lapse 
rate, the excess enthalpy was transported unward into the 
next higher layer. This process repeated itself until a 
zero- gradient of potential temperature was established in 
all layers involved in convective overturn. 

The National Space Science Date Center (G.S.F.C., 
Greenbelt, Md. ) provided the specific data sets of the 
IRIS, IR-Radiometry , and S-band occultation measurements. 
Magnetic tapes and microfilms were made available from 
which the particular dates were retrieved. The project 
utilized the rather unique occasion that the Mariner 9 ex- 
periments monitored Mars during three consecutive periods, 
i.e., (a) when the planet was enshrouded by dust, (b) in 
the transitional phase, and (c) after clearing. By 
running backward in time, the expanded computer model was 
attuned to Mariner 9 data first under dust-free conditions. 
Once this attunement was achieved . extrapolation into the 
transitional and fully dust-laden phases followed. 
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The basic o bje ctive of the investigation, reported 
herein, is to attune an expanded radiative -conductive - 
convective heat transfer mode l of the (CO 2 and dust) ab- 
sorbing and scattering Martian atmosphere to IRIS, IRR, 
and S-band occultation data , individually and concertedly, 
for the purpose of generating highly time -dependent at- 
mospheric and subsurface temperature soundings . These 
sets of tautochrones provided the basis for a reinterpre- 
tation of the aforementioned data as they permit to approxi- 
mate the vertical temperature distribution on Mars . 

A second objective is to contribute toward elucidating 
the interaction of Martian dust-aerosol and thermal energy 
transfer processes. 
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2.0 ANALYTICS OF RADIATIVE TRANSFER IN A MIE-SCATTERING 
AND ABSORBING MARTIAN ATMOSPHERE 
2.1 BASIC CONSIDERATIONS 

For 1 the transfer of radiative energy through a 
planetary atmosphere with vertically heterogeneous 
dust distribution , the physical processes of interest 
are the (particulate and molecular) absorption and 
multiple Mie -scattering of down- and upwardly directed 
solar radiative flux and the absorption, scattering 
and emission of Martian atmospheric and surface heat 
radiation under the condition of a given inhomogeneous 
stratification of aerosol and gases. 

In order to determine how the thermal structure on 
the planet Mars Is modified by a given aeolian par- 
ticulate load, one must compute the flux divergence 
for a specified number of layers, in atmosphere and soil, 
and a wide spectral range which Includes both short - 
and long-wave radiation. Although the solution of the 
radiative transfer equation for a specific atmosphere 
Is principally known (Chandrasekhar, 1960; Sekera, 1963a, 
b; Deirmendj ian , 196S) s the problem of computing -time 
economy arises. The exploration of the radiation field 
to any degree of completeness and accuracy would lead 
to an unacceptably large computational load (Braslau 
and Dave, 1973a, b) . 
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Recent.ly , Canosa and Penafiel (19 73) have shown that 
the spherical harmonics approximation to the radiative 
transfer equation represents a highly efficient method of 
numerically solving the transport problem for atmospheric 
models with highly anisotropic phase functions . The 
essential limitation in this approach is the restriction 
to homogeneous atmospheres. In 1974, Dave and Canosa dis- 
cussed a direct solution of the transfer equation for 
heterogeneous atmospheres by the spherical harmonics approxi- 
mation. They applied this method to atmospheric models with 
arbitrary vertical inhomogeneities due to variations in the 
scattering and absorption characteristics of the atmospheric 
constituents . In their formulation a the spherical harmonics 
equations represent a two -point boundary value problem for 
a system of first-order, linear differential equations with 
variable Coefficients . Through an extensive study , these 
analytics have been adapted to the needs of our project. 

2.2 GENERAL MATHEMATICAL FORMULATION 

In the two-point boundary value problem, Dave and 
Canosa (1974) introduced as independent variable the optical 
depth, measured from the top of the atmosphere, in the form 
of a linear combination of three distinct optical depths, i.e., 

(1) 
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The normal Rayleigh-scattering optical depth , T^ S ,r ' ) 

at any atmospheric level can be taken to be proportional 

to the pressure at that level. Values of this optical 

depth at the bottom of the atmosphere , ^ s ,r \ are given 

by Howard et alii .(1961) . Values of the Mie-scattering 

and particle-absorption optical depths, and 

respectively, can be obtained by following Deirmendjian 1 s 

procedure (1969). Details for computing values of the 

(a) 

CO 2 ” gas optical depth ,■ pertaining to the Martian at- 

mosphere, have been developed previously (Pallmann, 1968* 
Pallmann et alii, 1973). Obviously, the optical depth is 
highly wave-number dependent. 

The appropriate starting equation of radiative transfer 
for a general phase function. 


L 



is offered by Chandrasekhar (1960 , on page 150, eq. 92) with 
m = o (average intensity 1^°^) and making the coefficients 
in the Legendre polynomials series expansion in eq, (2) 
space-dependent. The average intensity is expanded in 
another finite series of Legendre polynomials with (L + 1) 
terms , i. e . ■ 
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and £jl representing the cosine of the zenith angle measured 
with respect to ■ the positive t“ axis (inward normall). By util- 
izing the orthogonality property of Legendre polynomials, i.e., 



one obtains an intermediate form of the transfer equation 
which if multiplied successively by the Legendre polynomials 
, wi 

yields the spherical harmonics equations. In matrix notion, 
they can be written as 


.th 


t- 


0 -, 1 = 


L, and integrated over 





<C f cr) 


<&<Lr) 


The appropriate, boundary conditions were developed by Marshak 
(19 47) and have the fo^m. 
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Dave and Canosa (1974) offer an explicit discussion of 
these expressions in the Appendix A of their paper. In 
addition, they explain in the main text of their paper 
what the bold-face letters in eq. (6) signify. For 
convenience , they are restated here * 
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0 <&/$ o dfs 
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and (It'D being the flux of a beam of solar radiation , per 
unit area normal to it . incident at the top of the plane- 
parallel atmosphere in the specified direction C C-ofo) 
as stated by Chandrasekhar (I960, p. 20, 145, 161, 216). 

The spherical harmonics equations (6) and (7) can be 
transformed by the finite difference method into an 
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equi valent algebraic problem which is particularly suited 
to treat transfer in nonhomogeneous atmospheres (see also 
Dave, 19 7 5a t b). It is worthwhile to mention that the 
adoption of vertical heterogeneity poses additional problems 
which are not taken care of by iust assuming the atmosphere 
to be divisible into a given number of layers each with con- 
stant radiative properties. In their aforementioned paper, 
Dave and Canosa emphasize that their algebraic problem is 
essentially different from that presented by Canosa and 
Penafiel (1973) although the method of solution is almost 
Identical. 


Considering some of the computational aspects, Dave 
and Canosa discuss the spatial structure of the model atmos- 
phere which is very similar to that adopted earlier in our 
research project. The upper- limit parameter L in the 
series of eqs. (2) and (3) is prescribed in practice through 
the numerical evaluation of the magnitude of in the 

(tr) matrix. If dust is absent 0^=1 for all £>2. Con- 
tinuity of the column vectors |p is assured by choosing L 
odd (or equal to 3 for dust-free conditions). 

By inspecting eq . (3) above, one realizes on the one 
hand that a very small number for L introduces strong os- 
cillations which are not realistic. On the other hand, a 


large numerical value for L will substantially enhance the 
computer time. Time calculations have been made which 
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permitted to optimize the problem. It is helpful to recog- 
nize that the spherical harmonics method is far less 
dissipative in nature than the iterative scheme for the 
transfer equation. 


2.3 SPECIFIC ANALYTICS OF A LOW-COST SIMULATION 

Chandrasekhar (I960; sec. 2.2) defined the downward and 
upward diffuse fluxes at a level , which we adjusted to 
the condition of the normal pointing downward, i . e . , 


Jo 


(9) 

Fgo =. 

and 

CIO) 

£«•>- 




The expressions indicate that these two flux functions which 
determine the spectral net flux at the level can be 

found by just computing the azimuth-independent This 

is an important simplification for any simulation procedure 
directed toward the determination of the net flux divergence 
and, finally, the time-rate of change of temperature at 
given levels 

A numerical value for the parameter L was adopted, i.e., 
L=ll. Although this might seem to be rather oscillatory, it 
is through the integration in the flux expressions (9) and 
CIO) that a smoothing process is applied. A systematic 
smoothing technique is available for various types of scat- 
tering (Kourganoff , 1952 ; Dave and Armstrong , 1974), which 
was not adopted in this project. 
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Sin ce the Dave- Canos a algorithm permits to treat 
heterogeneous conditions , the actual computation of the 
solution to the algebraic problem (Dave and Canosa's eq. 
(13) v 1974) is a much more involved problem. The numerical 
computation is unstable against the accumulation of 
round-off errors. The instability can be eliminated by 
carrying out stabilizing transformations on the matrix 
and vector sequences , involved in the algebraic system, 
at several conditioning noints. Its number was chosen to 

be N c ^ fi/d-H) + 3 t with L = 11 and o<t, <3 , the 
10 D 

number N of conditioning points ranged between 3 and 6. 
c 

The specific analytics for a low-cost simulation were 
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I 'cV - 1 fa 


a. >e 


jPszQ 

' I W ^ 

i 

as the azimuth-independent spectral radiance , 

M 


(12) 


F,c*r 


^ * o 


as the spectral downward flux, and 

ii 


(13) 


gcrj = 

Jj 


as the spectral upward flux. Marshak r s boundary conditions, 
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as discusser] by Dave and Canosa (1974; Appendix A) , can 
readily be specified for L » il • 

Computer time estimates indicate that for L - ii , the 
spherical harmonics approximation is one to two orders of 
magnitude faster on the equivalent of an IBM 360 computing 
system, model 195 with FORT RAN- H compiler, than the Gauss- 
Seidel iterative technique, A single precision generates 
numerical results accurate to the first 4 figures. In- 
corporating this algorithm into our RCCT-model and running 
the enhanced program to generate, over, a- period of 24 
Martian hours, 192 temperature soundings require several 
hours of computer time. 
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3.0 REFINEMENT OF TRANSMISSION FUNCTIONS IN "FLUX-OUT" 

AND "MARSTEMP" PROGRAMS 

Continued work with the C0 2 -absorption-line data, coin- 
piled by McClatchey et alii (1973), facilitated an essential 
improvement of the quasi-random OOj-transmission function in 
both the FLUX-OUT (F/0) and MARSTEMP (M/T) programs. It 
was found that band transmission modeling was superior to 
single- line approximation because of the well-known diffi- 
culties in line-overlapping, - asymmetry, and - shape factor 
(Kondratyev., 1965 , 1969 ; Goody, 1964), Consequently, an 
appropriate simulation of band absorption for Martian- 
atmospheric CO^ was computerized (sect. 3.1). An improved 
absorption capability was modeled which better represented 
intermediate line-intensities (sect. 3,2). Temperature- and 
pressure dependence of half-width , line Intensity, and shape 
factor found renewed attention (sect. 3.3) as well as the 
computational arrangement of spectral increments (sect. 3.4). 
Finally, the proper fitting of the band model to gaseous 
absorption was made utilizing the measurements of the 
Mariner 9 IRIS and IRR experiments (sect. 3.5). 

3.1 BAND ABSORPTION MODEL 

Goody (1964) has offered a rather detailed discussion on 
absorption band models. For computing numerical problems, 
it is desirable to develop averaging techniques which smooth 
over the finest structure in the bands without distorting the 
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band contour. Previously, we used the Ladenberg and Reiche 
(L-R) function as the basic mathematical form for CO 2 ab- 
sorption . Since we adopted a quasi-random model (Pallmann 
et alii, 1973), the L-R function appeared in the negative 
argument of an exponential function representing the trans- 
mission function. 

More recently, it was found that this type of trans- 
mission function does not optimally represent transfer of 
radiative energy under Martian atmospheric conditions since 
it does not take care of the many weak lines present. In 
order to correct for this shortcoming , another transmission 
function was adopted which made allowance for independent 
lines of unequal intensity and many weak lines before the 
strong- line region was reached* 

Basically, the analytical form of this transmission 
function is derived as follows : 


Assume a given wave-number range of N lines and each 
of these M lines as one member of an infinite array of 
randomly spaced lines. Let p($).dS be the probability that 
a line has an intensity between 3 and (S + dS ) . The adopted 


distribution function is 


(14) 




The cumulative probability of lines with intensity greater 


than S is 
( 15 ) 
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Laboratory measure merits have indicated that eq. (15) is 
fairly well satisfied. This indicates that the expression 
(14) accommodates a large number of weak lines (Godson, 1954). 
The statistically averaged absorption function reads 


AG*; An. | f(&) ~ *§•,) - 



for the Lorentz-line shape and (oO being half-width , 

“ S 3l ^ a _ 2^ 

u “ *. ■ ■■ ^ r— (a: amount of absorber mass), and /\n = 5 era 

A ‘fT’ot 

The designations I Q and 1^ stand for the zeroth and first 
order Bessel functions of the first kind with imaginary 


arguments'. 


The mean transmission function averaged over N super- 
posed arrays of lines is j ar **. 

(i7) T- 


A 


The restriction due to the Lorenz line shape is discussed 
in Sec. 3.3. 


3.2 MODELING INTERMEDIATE -LINE ABSORPTION 

The quasi-random transmission function, as given by 
eqs . (16) and (17), has two asymptotes, i.e., the weak- line 
and the strong-line approximations, for small and large 
u- values , respectively. Previously, we used this dual approach 
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to model "the whole u-range. A threshold value for u was 
prescribed, and the computer determined which aTvoroximation 
to use. Goody (1964* sects. 4.2.1) demonstrated bv means 
of the Matheson diapram that the intermediate u-range 
suffers fairly large inaccuracies. 

In applying this line of arguing to Martian atmos- 
pheric conditions, we decided that the intermediate range 
from 0.05 to 50.0, of optical path for the line centers 
needed further attention. In eq. (16), the modified Bessel 
functions with imaginary argument lent themselves to a 
direct numerical evaluation. Lewin and Grosberg (1952) offer 
a series expansion for this type of Bessel function, i.e., 


(18) 


with 


r 4 oo 

7 


i 1 j < 

f vJ C - <*•) — 

t 


GO 


a) 




^7 4' 


i: imaginary unit-J-1 

: gamma function 

Particularly, one obtains for the zeroth order function 

- .44 


(IS) 


I w-.TisL 
• fe 7?!r 


Klerer and Korn (1967) have given a numerical approxi- 
mation combined with estimates for the truncation error. 
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Computer runs were made for the purpose of determining the 

optimal form of representation of T and I. in the trans- 

o 1 

mission function. Acceptable forms were used first in 
our F/Q-program for dust- free conditions. The simulation 
of the IRIS spectra was quite satisfactory. It was in- 
ferred that the former combination of weak- /strong- line ab- 
sorption only was an Insufficient means of representing 
CO^ band absorption. The intermediate- line representation 
appeared to be a necessary refinement without which rather 
large differences between measurement and simulation of the 
IRIS spectra occurred. 

The subsequent incorporation of the Intermediate- line 
absorption capability of the interval 0.05£u450.0, into 
the M/T-program is considered to be one essential procedure 
in the attunement of the model to Mariner 9 data. 


3.3 TEMPERATURE AND PRESSURE DEPENDENCE OF LINE INTENSITY , 
HALF- WIDTH AND SHAPE FACTOR 

In McClatchey’s absorption data compilation, the line 
intensity S(T g ) is given for all known CO^ lines and T s =296K. 
Without the induced emission term, the temperature- depen dent 
line intensity S(T) at temperature T reads 


S(T) 


S(T s ) QUTJQrCt;) x 
QjiT) a r (T) 



( 20 ) 
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with 


-1 

E" : energy in cm of the lower state of transition, 

Q v , Q^: vibrational and rotational partition functions, 

respectively. 

Numerical information about the vibrational and rotational 
partition functions of COj (isotope 6-2-6) has been provided 
by McClatchey et alii (1973; Table 2). In addition, Goody 
(1964) has discussed some asnects surrounding the rotational 
partition function (in his sect. 3.2). On the basis of these 
discussions, it can be shown that eq. (20) may be replaced 
by 


( 21 ) 


SCO = stt) 



hmsE’Xr-Ts) 


r.r s 


This formula was subsequently used to compute integrated line 
intensities S for C0 2 at given temperatures T. Expression 
(20) indicates that S does not depend on pressure. 

The line half-width at half maximum, , depends on both 
temperature and pressure. KcClatchey et alii (1973) suggest 
the following relationship , 



with 

T = 296K; p = 1 atm; and (X : Lorentz line with at standard 

5 S 
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conditions. Although the temperature dependence T /T is 
not certain, measurements made with a CO^ laser (Ely and 
McCubbin, 1970) . indicate a linear factor (T /T) for the F20- 
line of the 10.4 CC >2 band. 

Concerning the line shape , some uncertainty exists. In 
the aforementioned absorption line compilation, the Lorentz 
shape was used. Due to Goody (1964) , the spectral absorp- 
tion coefficient can be expressed by 

(23) £0) = fu^»> = < 


where the Lorentz shape factor is normalized, i.e 


(24) 


(V (‘n-Hjd =. i 


-1 

with n: wave-number, cm „ 

The true line shape is not always given by the Lorentz 
shape which depends on pressure as discussed above. For 
Martian atmospheric pressures of less than 10 mb, the Doppler 
line shape is not negligible. The Doppler effect is associ- 
ated with molecular thermal motion. If an absorbing or 
emitting molecule has a velocity component w In the line 
of sight, a frequency shift A n - ^ .n Q about n Q is generated. 
As a consequence, the Doppler line shape adopts the form 


(25) 



with the Doppler line width being 
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(26) 



2hi 

c 



-i 

<jt c_m 


Calculations with numerical values , typical for Martian con- 
ditions , indicate that the Doppler width has about equal 
footing with the Lorentz width in the lower atmosphere 
gaining gradually with elevation. This is indicative of a 
combined Lorentz and Doppler broadening. Nevertheless, we 
have decided against the usage of such a combination since 
due to Goody (1964) the mixed Doppler and Lorentz line shape 
is too complex for analytical solution. Furthermore, the 
far wings of a line are always influenced by molecular 
collisions; if the absorbing path is long enough absorption 
will finally be controlled by the line-wing profile. Con- 
sequently the absorption will be given by eq. (16). 

Recently, it was pointed out by Rodgers (1976) that if 
both Doppler and collision broadening are significant the 
detailed line shape depends on the nature of the collisions. 
For most atmospheric gases, the Voigt line shape is adequate. 
Assuming that the two broadening processes are un.correlated, 
one averages the Doppler shifted Lorentz lines weighted with 
the distribution of velocities. As a result, one obtains a 
convolution of the Doppler and Lorentz shapes, i.e., 


(27) 



oC 



£ 

/ v4i ^.4 

-H- 4- n o ) + 
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For the numerical evaluation of the Voigt profile, see Kielkopf 
(1973). It is planned to accommodate this profile in the near 
future although no efficient numerical approximation exists. 

3.4 ARRANGEMENT OF SPECTRAL INCREMENTS 

The numerical treatment of the random transmission func- 
tion (see eq. (17)) requires each time an appropriate averaging 
procedure for N lines taken from the CO^ absorption bands. 

Since the project has been utilizing the McClatchey data com- 
pilation (1973), we had to accommodate about 9,900 lines in 

—1 

the far infrared between 200 and 2,000 cm and about 22,700 
lines in the near infrared between 2,000 and 10 ,000 cjn~^. 

The F/ 0-program which simulates the upwelling thermal 
flux at the top of the atmosphere, advantageously operated with 
a spectral increment of 5 cnr\ However the M/T program needed 
a more economical structure for these intervals because of 
the frequent iterations inherent in the time integration (of 
more than 144 time steps/day). After a careful inspection of 
the absorption line data, it was found that grouping of C0 2 

lines in intervals of variable length, generally between 20 

-1 

and SO cm was feasible if "lopsided" crowding of lines near 
either boundary was avoided. The new array consists of 140 
spectral increments appropriately covering the near and far 
infrared. 

Instead of the former spectral narameters, which were 
based on the Stull, t" T yatt and Plass (196 3) data and calculated 
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by Prabhakara and Hogan (1965), a new set of parameters were 
analytically developed and numerically determined . The data 
base for this set was found in the McClatchey compilation. 
From eq. (16) and the respective expressions for strong- and 
weak-line approximations, one obtains in accordance with 
Goody's procedure (1964) 



of lines, and d the average spacing . By inserting eqs* (21) 
and (22), we arrive at 


(30) XLzz 4 O. 


and 
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C 31) 


y 







line intensity at standard pressure p = 1 atm and 

vl 

temperature T = 296°K, in cm" /molecule, cm , as 
presented in McClatchey compilation ; 

-1 

Lorentz line ividth parameter at p , T , in cm /atm; 

s s 

and 

=1.4388.E“ energy of the lower molecular state in 

-1 

cm . 


■u 

The temperature factor (T =/T)^ represents the de- 
pendence of ot on T in eqs. (31) and (30). The additional 

ratio (T/T ) is derived from the temperature- dependent ex- 
s 

pression for the line intensity (see sect. 3.3; eq, (20) and 


( 22 )). 

The summations over the ilmdex "i” in eqs. (30) and (31) 


would represent a formidable task of computing if they are 
an integral part of the time-dependent M/T (or the F/0) 
program. Since these terms are exclusively temperature- 
dependent with the S^(i) and o^(i) to be read in from the 
McClatchey compilation, they were submitted to a polynomial 
approximation. Such a procedure admits a numerical corn- 
nutation, separate from the main program. Subsequently, the 
anproximated terms were stored on tape to be called in to 
the central processing unit whenever needed. 
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as 


C 32 ) 


and 


(33) 


For convenience, the terms in question are relabeled 




* 1 



where the symbols signify the same variables as before. 

These two functions were numerically evaluated for all wave- 

number increments which have the index "i" for computational 

convenience. (This was a re-indexing. ) Obviously, the 

values of F . and G . are not only different for each wave- 
0,1 0,1 

number interval, they are also temperature-dependent. 

The method of approximating by polynomials was deter- 
mined as follows: 

- Generate a second-order polynomial, dependent on tempera- 
ture , to approximate In M and In 

- Find for all Martian temperatures from 140°K to 

320°K and for each interval "i” , A^ln^F .land 

\ o,il 

BS= m (g C j .).. 

The data for all C0 2 -lines in a particular i- segment 
were read- in each time from the McClatchev tape. The pro- 
gram- run oroduced for M/T 137 segment values, and for F/0 
171 segments , of both ln(F Q ^ and ln^G Q ^). 
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Rub sequent ly , second-degree polynomials P(x) were fitted 

to both ln^F o and ln^G Q with the independent variable x 

in this case being 1.4388-CT ^ - T“^}. All numerical values 

s 

of these logarithmic functions were written in disk file for 
both the M/T and the F/0 programs. With the specific form 
for P(x) being 


(34) =? *- 

the functions F . and G„ . 

0,1 0,1 

' 35 > (Vi or (t,i 


4rX, +- C. X 4, 


are approximated by 



for all spectral increments of index ,T i ,T . 

In order to test the accuracy of the polynomial approxi- 
mation, the expressions 



were calculated and printed for all spectral increments and 
temperatures in the range 140K^ T^: 320K, with y being either 
ln|F . or ln[G The numerical values of the MAX- 

parameters average about 3%. The largest values usually 
occur at very low or very high temperatures which in turn 
appear rather infrequently in the Martian temperature pro- 
files. Therefore, it was concluded that the polynomials 

P(x) approximate the functions F . and G . with sufficient 

0,1 0,1 

accuracy. 
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3.5 FITTING BAND MODEL TO GASEOUS ABSORPTION CAPABILITY 

ALONG INHOMOGENEOUS PATH 

The. .analytical . expressions (30) and (31) of the 
previous section imply a coupling with the strong- artel weak- 
line forms of approximation. However, this procedure of 
anchoring the intermediate-line formula in these two domains 
requires a quantitative refinement by means of a numerical 
comparison of various band absorption models. Goodv (1964) 
has discussed the changeover from linear (u<^gl) to square- 
root (u>3) law regime by utilizing the M atheson diagram. 

From his discussion, it appears that a scale-shift in the 
variables u and y permits one to obtain a best fit. 

The numerical values of the shift factors were deter- 
mined through a process of attuning our model to the measured 
IRIS spectra of the dust- free period after January 1972. We 
found that a best fit was achieved by multiplying the u- 
parameter by (.25) and the y-parameter by (4.0). Various 
other trials of fitting did not generate as close an approxi- 
mation as that. 

The vertical heteoropeneity of the atmosphere required 
an additional refinement which was to he applied beyond the 
modeled stratification of 51 atmospheric layers inherent in 
the existing computer program (Pallmann et alii, 1973). Tt 
was decided to make use of the Curtis -Godson approximation 
as discussed by Goody (1964). The usual aim of scaling. 
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approximations in the case of non "homogeneous path is to 
reduce the problem to that of a homogeneous path with scaled 
amount a, scaled temperature T, and scaled pressure p. The 
scaling adjustments turned out to be minor since they were 
confined to the respective layer depths, usually 1 km or 
less (near the ground). Since we recognized the importance 
of lines of intermediate strength, we were left with an un- 
certainty in the scaled amount a of the order of (py^p)^ 
^ 1 . 02 . 

In conclusion, it was found that for the dust- free con- 
dition on Mars our model was closely attunabie to the measured 
IRIS spectra. The further step of attuning the model to the 
dust-laden situation will be discussed in chapter 6.0. 
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4.0 SIMULATION OF EFFECTIVE ABSORPTION BY MIE PARTICLES 
(DUST) 

In Chapter 2.0 of this report the analytics of a 
time-saving computer simulation was specified for the 
radiative transfer in a Mie-absorbing-scattering atmos- 
phere. However, a detailed time estimate indicated that 
the incorporation of the scalar MIe- algorithm Into our 
Radiative- Convective-Conductive Transfer (R.C.C.T. >. 

Model would need computing time for a Martian diurnal 
cycle of 10 hours on the CDC 6600 System 175. To achieve 
cyclic balance over 24 hours at least twice as many 10 
hours of computing time are needed. By repeating this pro- 
cedure for 3 different dates, the time requirement would 
rise to more than 60 hours. In the light of such a steep 
estimation, the project opted for a lower-cost simulation. 

The objective of the simulation is to concentrate on 
the impact of optically dense dust layers on the radi- 
ation field, and resultantly on the temperature field, 
rather than to formulate a sophisticated, but costly, 
model. Such a model would quantitatively describe the 
manifold Interactions of a heterogeneous CC^-gas and sili- 
cate dust stratification with solar and Martian-planetary 
radiation. The underlying analytics , as presented in 
Chapter 2.0, served as a framework for the approach chosen. 
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Various investigators have come up with findings per- 
tinent to the approach of elucidating the impact of dust on 
the aforementioned fields. Andres (1971) offered some in- 
sight into the mechanics of dust- lifting on Mars. He hypo- 
thesized that aeolian dust particles may have rounded shape 
because of perennial rubbing. Roundedness makes them amena- 
ble to Mie theory. Furthermore, Andres showed that the tur- 
bulent flow near the surface enhances dust lifting and sus- 
pension. Thus, thermal convection, as built into our model, 
may maintain a kind of randomness which justifies the choice 
of a Gaussian dust distribution function. The manner by 
which this Gaussian is related to the dust extinction co- 
efficient, as a function of height, was discussed earlier 
(Pallmann , 1974). 

Subsequent studies of dust saltation, erosion, and sus- 
pension by Sagan (1973) established findings consistent with 
those of Andres, but they failed to explain the great eleva- 
tion of several scale heights, at which dust was found during 
the gigantic storm event of 1971. Ajello et alii (1976) re- 
ported that limb observations by the Mariner 9 UV spectrom- 
eter television camera experiments showed that the dusty at- 
mosphere was optically thick up to an altitude of 50 km with 
a sharp cut-off and an extinction optical depth of 3, At 
times, the opaque dust layer was topped by an optically thin 
haze layer at 70 km. Pang and Hord (1973) estimated a mean 
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radius of about 2^m for the dust up to 50 km. Most recently. 
Pang et alii (1976) gave as mean radius a value of Ijajti. 

Goody and Gierasch (1973) created a model to simulate 
the incipient stage of a Martian dust storm. They proposed 
the local generation of free-atmospheric heat via radiative 
heating in an initial dust cylinder of critical size and sub- 
sequent synoptic- scale convective circulation of a vortex 
(dust hurricane). A dust devil of gigantic extent would be- 
gin to suppress turbulence favoring an organized vorticity 
field (Goody and . Gieras ch , 1974). The regional lifting of 
dust by such a circulatory system is likely to be an efficient 
means of generating high dust concentrations near the top of 
the atmospheric convection layer (tropopause) which may be 
more than one scale height in elevation. 

As shown by the hourly photography of the International 
Planetary Patrol (Martin, 1974) the dust storm started locally 
in the NE periphery of Noachis with another center in the NW 
quadrant of Hellas. A third dust -^mushroom 11 appeared at about 
the same latitude over Ausonia but stayed only intermittently. 
The photographic evidence supports Goody and Gierasch 1 s model. 
If one assumes the existence of powerful jet streams in the 
vicinity of the tropopause, the dust would spread enshrouding 
the globe as fast as it did (15 days). 

The combination of gravitational settling and vertical 
mixing (Conrath, 1974; 1975) contributes toward explaining the 
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observed time evolution of the silicate particulates in the 
atmosphere. It is, however, not necessary to assume intense 
vertical mixing in all atmospheric layers throughout the dust 
stratum and for all longitudes and latitudes. As discussed 
in Chapter 7, near-adiabatic lapse rates which appear as a 
consequence of vigorous convective mixing, may be preferen- 
tially present in the lowest scale height and, again, near the 
top of the dust. In the center of the dust stratum, say, 
between 15 and 30 km, the atmosphere might be rather stable. 


4.1 SOME ANALYTICS AS BASIS FOR A LOWER-COST SIMULATION 
OF THE DUST IMPACT 

In equation (2) of Section 2.2 of this report, the normal- 
ized phase function L§_. 

(2) 


C_0-3 o) — 'o'j <£w o) 

o 


is presented as a series expansion in Legendre polynomials 
$^(cos(0) with (&) being the scattering angle. The series 
coefficients may be represented by 

( 37 ) 'ey = — — 

6 ir/3 


( $<, 


iresenrea ay 
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J_. L^.W + l— (x,^)yi(t)dr 
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In this form, the 


CCanosa, J, and H. R. Penafiel, 1973) 

following designations are used: 

t : radius of aerosol particle, 

3 

'YlCr): number of particles per cm per unit radius interval at the 

radius Y“ 


X; 


ofrfTr 

X 


y Mi e- factor, 
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fj : volume scattering coefficient. 

< i ) 

The Legendre coefficients have been treated in detail 

by Dave (1970). The expression (37) requires explicit knowledge 
of the size distribution function n(r) and the scattering 


coefficient 


Psc' 


As discussed in Section 4.3, we showed that 
Deirmendjian 1 s (1939) Silicate Haze M model is best suited to 
represent a dust load consistent with the quantitative in- 
formation on vertically integrated column number density N, 
and mass density m as offered by Conrath et alii (1973; Sect. 
E, Table XX-4). The particulate size distribution function is 
thereby modelled as 

(38) x(T) ~ 3333 K 

3 

with the total number per cm of particulates being N - 100. 

Due to Deirmendjian* the volume scattering coefficient 
may be written as 




A 


00 


(39) 


pJA, ncr)) = trj r * ncx) [XJi 


pC c rj J,: 


where K designates the scattering efficiency factor C called 
s c 

Q sc by van de Hulst, 1957). Incidentally* the integral on the 
right-hand side creates the bias (r ) in favor of the larger 
sizes in the spectrum. For Mie-factors x-^5, K (r) under 

5 C 

the integral may be replaced by K 2 . Thus, one obtains 

s c 


(40) 

with 

(41) 


jSgtCA/KCT*)) =* Yi,C~r) 


d/r 


f ncr)d 





Since (r ) is a non- decreasing function of r, we may write 
(42) -n(y) <Lt =S (rf) Jhct) c(r = (r*). fsj 


Finally, this procedure renders 

( 43 ) = £tr(r 3 -)N 

This relationship will he utilised later on. 

In order to improve on our first-order dust transmission 
function which we formerly attached by multiplication to the 
gaseous transmission function, we developed an analytical ex- 
pression taking up the additional absorption of radiative 
energy via multiple scattering by dust* Tbe underlying concept 
is that multiple scattering by dust enhances the radiative path 
length through the atmosphere increasing thereby the amount of 
radiative energy refracted into the particulates and, subse- 
quently, absorbed* The Imaginary portion of the complex index 
of refraction regulates the portion being absorbed. 

We start with the scalar equation of radiative transfer 
within the spectral increment 


(44) 


n — T<X 

' dl'CV ■ * 


The notation is standard; for convenience, see Section 2,2 of 

in 

i 


this report. The source function 

.'tv 




(45) 


is composed of two terms* The first term on the right-hand 




-34- 


side 'represents the portion of the solar beam, in the direction 
attenuated over the slant optical range (£!^/ 
and subsequently scattered into the direction (^,-*jb). The 
second term represents scattered "self-illumination”, that is 
to say, the fraction of radiative quanta impinching as diffuse 
radiation from all directions ( jvl , ) which becomes scattered 

into the direction C The function ^(T^) is equal to 

($&>/&&>)- 

With the prevalent size of the Martian dust particle being 
between l^tm (Pang et alii, 1976) and 2 (Pang and Hord, 1973), 
a rather pronounced forward scatter is present* Each layer 
passes on to the next lower layer through this maximum of 
scatter most of the scattered energy* Given a manifold layer 
structure of dust load and multiple scattering therein, the 
diffused radiation field gradually approaches axial symmetry 
such that the diffuse intensity near the local zenith is just 
a few times larger than that near the horizon (see Sobolev, 1963; 
his chapt. 3). If we approximate this characteristic by the 
condition of isotropy, based on the principle of energy con- 
servation, we are able to specify the analytical form for the 
dust transmission function. 

The radiative transfer equation for 1^°^ is multiplied 
by dw and integrated over all solid angles. We obtain for an 
isotropic source function (Goody, 1964; p. 52: Eddington's 
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approximation ) 

fa cU v 7 

x: geometric depth* 

For computational purposes , we adopt here a two-stream model, 
A split of the flux function into F'^ and F^" appears appli- 
cable for a scattering and absorbing source function because 
of the initially pronounced forward scatter of Mie theory. 

We find . 


it 


CH7> 


4ir j3dx V/ ^° + 


©Alfdo/l 

(see also Deirmendjian, 1969; p, 94). 'fowuiarJ ~*pjt I 


For isotropic conditions , F^ = 7T , and p = const. 
Thus, one obtains 

UF* . 


(48) 


5* 


ijfitrf (Q/lijfLe. 




On the right-hand side, the first term represents the hemi- 
spherical downward flux extinction, the second and third 
terms the attenuation due to primary and multiple scattering, 
respectively. 

Integrating over a finite depth range /\x and inverting 
the natural logarithm yields 

i /?>■*& 


(49) 
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■= Sp” = -x^p **■*!>■[ 
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The second exponential factor on the right-hand side stands 
for the component transmission function due to primary and 
higher order scattering. Since the first exponential factor 
previously represented the only dust absorption function used 
in our model (Pallmann, 1974), the second factor has been 
taken up to account for the additional absorption of radia- 
tive energy due to multiple scattering. The replacement of 
the scattering coefficient |3 sc by an effective dust ab- 
sorption coefficient will be dealt with in sections 4.2 and 
4.3 of this report, 

4.2 BUDGETING IN BULK OF ENERGY ABSORPTION IN GROUND- 

ATMOSPHERE SYSTEM 

With a model silicate haze M (Section 4,1) and an albedo 
of single scattering ■'Sf = .8 for visible solar radiation 
(Deirmendjian, 1969$ Table T.99), we find a 20% absorption 
in the dust extinction process. In the top 10 km of the 
model atmosphere, from 50 to 40 km of height, there is no 
dust. The dust spreads out from 40 km down to the surface 
in accordance with a Gaussian distribution centered at ground 
level. Because of the proportionality between the number 
density N of dust and the effective dust absorption coefficient, 
a plot of the latter function versus height may represent the 
shape of the distribution curve. 

A tabular account may elucidate the budgeting procedure: 
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50-40 km 

No dust 


40-0 km 

Dust layer 


40-35 km 

Reflectivity of down- 

Transmissivity of down- 


ward solar radiation: 

ward solar radiation: 


<♦05) per 1 km 
total fraction 

( . 95) per 1 km 


reflected: (.23) 

transmitted: (.95) S = .77 

35-0 km 

Per 1-km layer: 



(.9934) transmitted 

(.0066) absorbed 


or finally at ground, 

in 35 km layer: 


(.9934) 35 = (.793) of 

(.77) - (.6107) - ,1593 


impinchment (.77) 
transmitted : .6107 

absorbed 

For solar 

radiation imp inching at the 

ground: 


Surface Albedo: ( 

.3) 

0 km 

(.6107) x C. 3) = .1832 

(.6107) x (.7) = .4275 


reflected $ 

absorbed by ground j 

0-35 km 

(.9334) 35 = (.793) of 

absorbed by atmosphere: 


ground-reflected portion 
( ; 18 32 ) transmitted: 
(.1453) 

(.1832) - (.1453) = (.0379) 

40 km 

Transmitted outward: 
(.1453) 



Total fraction of solar radiation absorbed in the down- and 
up-sweep : 
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.1593 
* .0379 
.1972 

or close to 20%. This is a satisfactory match to the value 

of Cl - CD' ) of the model silicate haze M whose complex index 
o 

of refraction is m = 1.55 - i* 0.0155. The planetary albedo 
amounts to < * 375 3 > which is rather high. 

A second budget with a planetary albedo value of (.29) 
yields the following fractions of reflection, transmission and 
absorption: 

5 0-4-0 km No dust 
4-0-0 km Dust layer 

4-0-35 km Reflectivity of down- Transmissivity of down- 
ward solar radiation: ward solar radiation: 

(.05) per 1 km (.95) per 1 km 

total fraction 

reflected: (.23) transmitted: (.77) 

35-0 km Fer 1-km layer 

(. 98094-) transmitted (.01906) absorbed 

or finally at ground: in 35 km. layer 

(.98Q9H) 35 = (.5099) (.77) - (.3926) = C.3774) 

of impinchment (.77) 
transmitted: (.3926) 


absorbed 
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For solar radiation impinching at the ground: 

Surface Albedo: (.3) 

0 km (.3926) x (.3) = (.1178) (.3926) x (.7) = (.2748) 

reflected y absorbed by ground ^ 

35 

0-35 km (.98094) = (.5099) of absorbed by atmosphere; 

ground-reflected portion (.1178) - (.0601) * (.0577) 
(.1178) is transmitted; 

(.0601) 

40 km Transmitted outward: 

(.0601) 

This time, the total fraction of solar radiation obs orbed 
in the down- and up-sweep through the atmosphere is 

(.3774) 

+ (.0577) 

(.4351) 

or close to about 44%. 

In summarizing the two bulk budgets of energy absorp- 
tion in the ground-atmosphere, we have in percent of solar 
radiation impinching at the top of the atmosphere; 
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Case with 

Transmissivity 

between 35 and Diffuse 


0 km being 

Reflectivity 

Absorpti vity 

1- 

37.5% 




Outer Space 

2— 

29. % 


1- ( . 99 340) /Jem 

23.% + 14.5% 

19.8% 



Atmosphere 

l-(. 98094) /km 

23.% + 6.% 

43.7% 

1- 

30.% 

42.7% 



Ground 

2- 

30.% 

27.3% 

Balance : 



1- (37.5 + 19. 

8 + 42.7)% = 100% 


2- (29.0 + 43, 

7 + 27. 3 ) % * 100% 


4.3 EFFECTIVE 

ABSORPTION BY MARTIAN DUST PARTICLES 


Although the value for atmospheric absorption given in 

the previous section seems rather high, it is consistent with 

the extinction optical depth of 3, offered by Ajello et alii 

(1976). Furthermore, the vertically integrated column number 

7 

density of 1.21 x 10 for r - 1 as given by Conrath et 
alii (1973), yields, when inserted in equation (43) 

pj x > N ) - 

a value for the volume scattering coefficient Q of about 

( SC 

1.2 x 10“ cm” . In this context, one has to remember that 

— ? — y v f 2 

(r ) must satisfy the relationship (r ) = jj**j r;n(r)dr. 





-41- 


This expresses the characteristic of a poly disperse dust 
particle size spectrum (Esirmendjian , 1969). 

The previously adopted principle that, via multiple 
scattering within an optically-thick dust layer, the likelihood 
increases for photons to become absorbed, mostly by the par- 
ticulates, has led to a numerical estimation of the percentage 
involved in this additional absorption. Formerly, we had just 
accounted for an average slant solar beam absorption in the 
mode of single scattering (Pallmann et alii, 1973; Pallmann, 
1974; Dannevik and Pallmann, 1974). Based on calculations 
of the brightness of the sky in the case of multiple scattering 
(Atroshenko et alii, 1963), it was determined that multiple 
scattering takes up about 25% of the total scattering leaving 
75% of the diffuse sky radiation at the surface for primary 
scattering. With an albedo of single scattering value of 
(.63), the ratio of the coefficient of multiple scattering to 
the extinction coefficient is close to (.21). 

Deirmendjian (1969) obtained for model silicate haze M 
and various values of the imaginary part of the complex index 
of refraction the rather even numerical figure of about 
1. x ItT 6 cm*" 1 for the volume extinction coefficient. With 
this particular value and adopting an approximate 50% effi- 
ciency for multiple scattering to generate further absorption, 
one arrives at (.5) x (. 21) «f(,10). Together with the rela- 
tive dust absorption coefficient, previously incorporated into 
the model, we have 



C5D) 




direct) 


or with 


p 




Jb (mult. so.) = ( ' ZB) 


ex 


t -, n — 6 -1 

= 1. x 10 cm 


p 


ex 


(51) 


P 


Dust 


ab 


= .2.6 x 10 


“7 


cm 


**1 


This value of the effective absorption coefficient due 
to dust should be considered as representative of atmospheric 
conditions close to the maximum dust-load in space and time. 
However , it has been derived from a vertical imp in chment of 
solar radiation rather than time-dependent slant incidence. 
The diurnal cycle of the apparent march of the sun in the 
local sky covers a zenith angle range from ('ft' / 2) to noon 
zenithal distance, equal to the difference between latitude 
and solar declination. If one takes in the average a zenith 
angle of (0.3f^ the corresponding value of sec (0.3frO for 
the slant range is close to ( 1.6). Therefore, a somewhat 
enhanced value for the maximum of the effective absorption 
coefficient may be allowed. Also, the gradual dissipation of 
the dust phase through November and December 1971 calls for a 
persistent decline of the absorption of radiative energy by 


dust. 


The circumstances just described suggested a series of 
simulations for the effective dust absorption coefficient EP 
in order to determine the numerical values of A and B in the 
expression 
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(52) E p( Z) =* A e*xp 3 Z j* 

In an earlier report (Pallmann, 1974) the compatibility 

of this analytical form with other forms (see Conrath, 19 7 4) 

was discussed* A series of computer runs with our R.C.C.T. 

and F/0 models was undertaken* The numerical values for the 

amplitude parameter A ranged from (0.5 x 10 ) cm to 

(3.0 x 10" 7 ) cm"'*'. The parameter B was held at (1.11 x 10 ^ ) 
-2 

cm .In chapter 7 of this report $ several CALCOMP plots of 
the simulation are presented. 
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5.0 IRIS y IRR AND S-BAND OCCULTATION DATA 

In order to perforin a multiple attunement of our 
time-dependent radiative-convective-conductive heat 
transfer (R.C.C.T.) model to various types of Mariner 9 
measurements, IRISyIRR-, and S-band occultation-data 
were requested from the National Space Science Data 
Center (N.S.S.D.C.) at Goddard-Space Flight Center in 
Greenbelt, Maryland 20771, 

The Data Center made two sets of magnetic tape 
available which held, the measured information of IRIS and 
IRR for several hundred spacecraft revolutions. In 
addition, the N, S. S.D. C. sent a set of microfilms of 
tabulated and plotted S-band occultation data for both 
the standard and extended missions. 

Because of the extensive volume of data on the tapes, 
a rather time-consuming computer sorting became necessary. 
Some problems with the retrievability of IRIS data blocks 
had to be overcome before the excerpts of the desired 
data were put on about 20,000 cards. We found that the 
actual IRIS measurements had been arranged in an excellent 
fashion. The N.S.S.D.C, Information Packet had been care- 
fully prepared serving efficiently its purpose of explaining 
the code of the -IRIS Reduced Data Record, Type 7- — Cali- 
brated Martian Spectrum. 

Similarly, the IRR tape data had been very well 
arranged. The aforementioned problem did not exist here 



-45- 


since the data blocks were labeled such that excerpting was 
straight forward. The microfilms with the S-band occultation 
data were handled with standard reader-printer equipment 
available in the Department of Earth and Atmospheric Sciences 
at Saint Louis University. 

A brief word on cost concerning the tape data sorting ar.d 
excerpting procedure. Because of the unlabeled formats used 
for the IRIS sequences, about $250— at University rate — • 
had to be expanded on an IBM 360/65. The cost of the IR- 
radiometric data handling was about $170 — whereas microfilm 
printing cost was minor. 

Furthermore, the attention of the reader is directed 
toward a peculiarity in the data, which is related to the 
absence of coincident observations among S-band occultation, 
IRIS, and IRR measurements. Orbital conditions of the space- 
craft caused, for a given region and local time, the IRIS ob- 
servations to lag behind the occultation measurements by 15 
to 45 days (Leovy, 1973). However, S~band occultation data 
can be used in conjunction with IRIS or IRR measurements if 
larger-scale Martian characteristics within limited latitude 
belts are of prime interest. 

An additional complementarity may be found in that 
occultation measurements, although noisier, generally possess 
a better vertical resolution than the IRIS information which 
offers reliable results over deeper layers. In fact, Conrath 
(19 74b) determined the vertical resolution of the technique 
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for the IRIS- inverted Martian temperature profiles. He 
offered a value of about 6.5 km in the lowest 10 km increas- 
ing gradually to slightly more than 20 km above the height 
of 33 km. 

The data card bank, g rated from the magnetic tape 
and film material, served as the base for an extended simu- 
lation procedure. It was finally decided to concentrate on 
the following dates for the purpose of representing the 
gradual transition In the S.H. mid- latitude from the dust- 
laden to the dust-free condition: 

— • November . 22 , 1971 (revolution 17). 

— December 03, 1971 (revolution 39), 

• — - December 24, 1971 (revolution 79), 

- — January 08, 1972 (revolution 111), 

— • February 16, 1972 (revolution 189). 

For the last 2 dates, no S-band occultation data were avail- 
able. They reappear after May 6, 1972 when the tracking 
station at Goldstone, California received occultation signals 
during the even-numbered revolutions from 352 to 450. 
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6.0 REFINED MODEL AND ITS ATTUNEMENT TO DUST-LADEN AND 
DUST-FREE CONDITIONS 

In the present chapter, the governing mathematical 
equations of the model are presented in section 6.1 as 
they have been refined on the basis of our former radia- 
tive-conductive model* A specific amplification of this 
model is given by the incorporation of convective heat 
transfer anywhere in the atmospheric stratification, not 
just near the ground* The principal refinement, however, 
should be seen in improving the CO^ and dust transmission 
functions. The former makes allowance now for inter- 
mediate-line strengths as well as strong- and weak- line 
absorption whereas the latter Includes indirect absorp- 
tion by optically thick dust via multiple scattering, 
adding to direct dust extinction. Some procedural guide- 
lines of the overall attunement of the model follow in 
section 6.2. This discussion deals also with the tuning 
problem for dust-laden conditions, 

6.1 GOVERNING ANALYTICS OF THE MODEL 

The solution of the radiative transfer equation for 
the planetary problem leads, when Integrated over the 
full range of (4'?r) steradian, to an expression of the 
spectral radiative net flux (Pallmann et alii, 1973). If 
additionally integrated over all wave numbers, the total 
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net flux results which may be inserted in the atmospheric 
heating rate equation to render 



In this equation, the subscript n m" designates the In- 
dex for a given level between the top of the atmosphere 
(at 50 km) and the surface, and the subscript U I" indicates 
a particular radiative wave-number increment. Furthermore, 
S^tO.) is the solar spectral irradiance at the top of the 
atmosphere, Cr^(x s ) is the upwardly directed spectral flux, 
emitted by the ground, at depth x , (B * is the Planck's 
function, $?. is the flux transmission function for CO-, 

X A 

is the dust transmission function (grey absorption), and k 
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is the thermal conductivity. The derivatives in this equa- 
tion are replaced by finite differences. A simple forward 
difference is used for the time derivative, while modified 
centered differences are employed in the space derivatives. 
The resulting equation is quite lengthly and will not be 
reproduced here. 

The sub -surface heating-rate equation in its finite 


difference form is given by 

■ -■ Q 


( 53 ) ^C. 


fs*. 


At 


A** 



f 


X 




or* Vlt ^ 


for qu = N *f~ ^ ; H t ; ■ - * ^ 


ct 


where the superscript tt n" refers to the time step, and 
is the lowest soil level under consideration. 


The ground- atmosphere interface condition at x = x, 

" J s 

which expresses the continuity of the total heat flux across 

thp i ntp“pfaGf> . is ar 








-50- 


or 7J(v-°; t; = T?v ^ ^ . 

In equation (54), 3c designates the soil thermal conductivity, 

s 

£ the bulk surface emissivity, <5* the Stefan- Boltzmann con- 
stant, and T g the soil temperature. 

In equation (53), as well as equation (54) and the finite 
difference form of equation (52), the DuFort-Frankel two- 
level explicit differencing method has been employed (Carnahan, 
Luther and Wilkes, 1969). This procedure is computationally un 
conditionally stable. 

The upper boundary condition for the ground- atmosphere 
system reads 


(55) 


T(o,-t] * T 0 


for t ^ O 


The lower boundary condition is 

(SB) TsC^'t) = 7j ) fat- t» 0 . 

The initial conditions are given by 
(57) ~ J f ofXs< * <Xd 

T«,oJ = T f C<J J for 0$x$x s 

as a discrete input * 

The system of equations (52) - (57) suffices to estab- 


lish a unique solution at 81 grid points in the vertical 
throughout the region ^ x < 0, for any time t> 0 
(30 grid point in the soil, and 10 in the lowest 6.5 km of 
the atmosphere). 
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The built-in conve ctive heat transfer capability is 
based on Kraichnan T s turbulent thermal convection theory 
(1982). Thereby, the convective heat flux at the top of 
the constant flux layer (CEL) is given < 

( 58 ) 
with 

c = 8.2 x 10 6 ergo k” 1 

J a 2-1 

K - 8 . cm sec 

, in 2-1 

V = 10. cm sec 

g = 3.76 x 10^ cm sec" ^ 

?CFL = !*i 2 x g cm" 3 



The mathematical forms of the CO^ and dust transmission 
functions have been presented in chapters 3 and 4 respec- 
tively. Attention is called to the fact that each of the 
transmission functions ( ) in equations (52) and (54) 

is approximated by a product which may be written as 


(59) 





such that the product index M j M extends over all J interven- 
ing segments (x^ - x j-l^ ^ e-fcween the end points of the 
depth interval (x^ , x^}. For convenience, these segments 
taken so as to coincidence with the atmospheric layers 
labeled by the index ’'m 1 ' . 


are 
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Questions as to the size of space and time increments 
as well as remarks concerning the previous model have been 
treated elsewhere (Pallmann et alii, 1973). 

6.2 HIGHLIGHTS OF ATTUNEMENT 

The procedure of fitting our refined band model to the 
CO^ absorption capabilities was covered in section 3,5 for 
an inhomogeneous radiative path. As indicated earlier, one 
of the goals of the simulation was to achieve a fitting as 
close as possible for the dust-free condition. The particu- 
lar date chosen to represent this condition was February 16, 
1972* IRIS and IRR measurements for the midlatitudes of the 
Southern Hemisphere had to be compared and adjusted on the 
basis of their respective resolution. Our R.C.C.T. model 
had to generate a sequence of simulation runs which produced 
temperature- soundings cyclically balanced over a period of 
24 hours (diurnal cyclic balance). Between 288 and 432 
time-steps of Integration or 2 to 3 days ware needed. A 
cyclically balanced temperature sounding which coincided with 
the time of the IRIS measurement, served as input for F/0- 
model rims whose output was a simulated IRIS spectrum over 
the wave-number range from 200 to 2,000 cm - ^. Depending on 
how close a match was obtained, the numerical values of 
various scaling parameters had to be adjusted. We refer 
here to the K parameter in equation (14) , the modified u- 
and y-parameters and the Curtis-Godson scaling factors of 
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section 3*5. Additional nuns of both the R.C.C.T. and F/0 
programs with adjusted, values for these parameters were 
undertaken. This iterative mode of attunement finally led 
us to accomplish a close match to the measured IRIS spectrum 
produced under dust- free conditions over the midlatitudes 
of the Southern Hemisphere. 

In turning next to the dust storm dissipation phase we 
selected the IRIS, IRR, and S-band occultatlon measurements 
of December 24, 1371 as input data for another series of 
runs with both programs. This time, a dust transmission 
function had to be inserted to model a moderate dust impact 
on the solar and planetary radiation fields. By carefully 
evaluating the outputs of time- dependent T-soundings, 
temperature lapse rates, and convective overturn from the 
R.C.C.T. model, and upward spectral flux, emitted by ground 
and atmosphere and transmitted to the top C 5 0 km), from 
the F/O-program, we adjusted the value for the effective 
dust absorption coefficient, atmospheric layer absorptivities , 
and cumulative reflectivity at the top. For numerical detail, 
reference is made to chapter 4. 

For a treatment of the optically-thick dustload in the 
Martian atmosphere, the day November 22, 1971 was chosen. 
Again, IRIS, IRR and S-band occultation data were scrutinized 
to be used as input for our simulation procedure. It became 
apparent that the amplitude value of the effective dust 
absorption coefficient had to be quite large. Vie adopted 
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-7 -1 

values as high as (3,0 x 10 ) cm which corresponds to an 

i 

absorption optical depth of Cl. 2) or an extinction optical 
depth of about (3.0). More detail may be found in chapter 4. 
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7.0 DISCUSSION OF RESULTS 

With the attuned computer model, as discussed in 
chapter 6 , a series of outputs was generated via the 
time-dependent MARSTEMP (M/T) and adjunct FLUXOUT CF/0) 
programs. The M/T program produced, for each time-step 
(between 4 and 10 min) over a total of 36 to 60 hours 
of simulated time, the following outputs as printouts: 

— IR-flux emitted by the ground and transmitted 
to 52 explicit atmospheric reference levels; 

— solar flux, downwardly directed, and transmitted 
to these 52 levels; 

— solar flux, diffusely reflected by the ground, 
and transmitted upward to these 52 levels; 

IR-flux, emitted upwar d by 51 atmospheric source 

layers , and transmitted to 51 reference levels; 

— IR-flux, emitted downward by 51 atmospheric 
source layers, and transmitted to 51 reference 
levels. 

— ■ upward total flux at 52 levels; 

— ' downward total flux at 52 levels; 

— net flux at 52 levels; 

- — new temperature at 52 levels; 

— rate of temperature change at 52 levels; 

— — temperature lapse rate at 52 levels; and 

— subsurface soil temperature at 30 depths. 
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The F /0 program generated for a given time (generally 
identical to the indicated MARINER 9 measurement hour) the 

spectral flux value within 360 increments of 5 cm -1 between 

-1 

200 and 2,000 cm (IRIS range). These fluxes were avail- 
able as 

— IR-flux emitted upward by the ground, and trans- 
mitted to the top of the atmosphere 5 
— IR-flux emitted upward by an individual atmos- 
pheric source layer, subsequently attenuated by 
the intervening layers , and then Integrated over 
all source layers; and 

— sum of the previous two fluxes for all spectral 
increments . 

From this broad variety of outputs, mostly the temperature 
soundings of the MARSTEM simulation and the F/ 0-outputs 
will be discussed. In addition, some plotted input functions 
will be presented as they relate to the findings . 

7,1 ATTUNED MODEL OUTPUTS, DUST FREE CASE 

The first basic step of attuning the model to Mariner 9 
data was to fit the refined gaseous CO^ band transmission 
function (chapter 3) to IRIS measurements taken under the 
condition of a dust- free Martian atmosphere. By means of an 
iteration of gradually Improving computative F /0 attunements, 
a rather close fit of the upwelling flux output to the 
measured IRIS spectrum was achieved. For February 16, 1372, 



FLUX 



WRVE-NUMBER , C 10/CM 

Fig, 1: Spectral upward flux, at the top of the Martian atmosphere, in erg, sec”i cm“2 

(hemisph. )~’Vcirf'i, versus wave number in cm-1 for February 16 , 1972 , 12.72 LMST, 
at -41. 4 deg. lat. and 327,7 deg, long., as simulated by unabridged F/0 program 
on the CDC 6 600 System 175 (marked by ,, + lt sign). Flux is the sum of ground and 
atmospheric layer emission, duly attenuated, — Solid curve: average of 3 

measured IRIS spectra in degraded mode. 
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the quality of the attunement is illustrated in figure 1. 

The matching upward flux was computed by an unabridged F/0 
program in which the series expansion for the modified 
Bessel function of zeroth and first order held 99 terms 
(see equations (16) (17), (18) and (19)). As mentioned 
earlier, the representation of CO^ band absorption for 
intermediate line strength became necessary, since our simu- 
lation study demonstrated that neither a strong-line model 
alone, nor strong- line-weak- line models rendered sufficiently 
accurate numerical attunement. It is exactly the attune- 
ment procedure which provides a sufficiency criterion. 

The temperature sounding used as an input for the F/0 
run of figure 1 was a cyclically balanced T-sounding gener- 
ated by the unabridged M/T program on the CDC 66Q 0/System 
175. (This computing system of United Computing Systems, 

Inc. , St. Louis is about 13 times faster than the CDC 3300 
at Yalem Computer Center, Saint Louis University, as far 
as our particular program is concerned. Because of limited 
computing funds, the simulation study with the system 175 
was kept to 4,5 hours. Additional computing was performed 
through the easily accessible Yalem Center at Saint Louis 
University). 

In order to throw more light on the question of how 
the number of terms in the aforementioned series expansion of 
the Bessel functions affected the output, some additional 
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Fig. 2: Same as Fig. 1 except that this output was simu- 
lated by the abridged form (33 terms) of the 
- F/ 0 program on the CDC 3300/YALFM; the measured 
data are an individual IRIS spectrum in a non- 
degraded mode. This time, the curve of the 
measured spectrum is marked by an asterisk; the 
simulation is given a solid line. 
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runs with just 33 terms were carried out. With this par- 
ticular modification, we combined some variations in the 
temperature input. Thus, instead of feeding into the F/0 
run a T-sounding which was generated by a run of 48 hours 
of the M/T model (cyclical balance), we used the original 
IRIS inverted temperature sounding in our F/0 simulation. 

Figure 2 shows the result. 

It is apparent that the matching is of good quality. 

Only minor differences exist due to the adoption of the 
abridged form for the series expansion of the modified Bessel 
function. The surface temperature and the T-sounding for 
the simulation was the same as that of the measured IRIS 
spectrum as shown in figure 3. Attention is called to the 
steep drop of temperature near the surface. Otherwise, 
there is nowhere a lapse rate close to the adiabatic value 
( 5 , 2K/km) . 

As an additional variation, the IRIS-inverted T-sounding, 
with the mean IRR brightness temperature of 259. 5K for the 
surface, was inserted in the abridged F/0 program. The out- 
put is presented in figure 4 which clearly shows the in- 
fluence of the surface temperature on the spectral flux welling 
up at the top of the atmo sphere , The respective input sound- 
ing is shown in figure 5. 

From the illustrations, previously shown, it is evident 
the surface temperature is a controlling factor in any 
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Fig. 3: The temperature sounding used as input for the 

simulation whose output appears in Fig* 2. 




HAVE-NUMBER, tlO/CMJ 


Fig. 4: Spectral upward flux, at the top of the Martian 

atmosphere, in erg sec~l cm - 2 Chemisph)*“l/cm“l, 
versus wave-number in cm -1 for February 16, 1972, 
12.72 LMST at -41.4 deg. lat. and 327.7 deg. 
long., as simulated (solid curve) by the abridged 
F/0 program with T-sounding of Fig. 5 as input 
(T = 259. 5K) * Curve, with (*) is the measured 
IRIS spectrum. 



Fig. 5: Temperature sounding used as input for simulation 

whose output is given in Fig. 4. Curve is identi- 
cal to T-sounding of Fig. 3 except for T s = 259. 5K. 
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Fig. 6: M/T-model generated temperature sounding versus 

height used as input for a F/ O-simulation under 
dust-free conditions on February 16, 1972; 

12*72 LMST. Areographic location; -41.4 deg. Xat. 
and 327.7 deg. long.; temperature lapse rate up 
to 6 km: adiabatic or steeper. 
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simulation procedure. It is for this reason that we have de- 
veloped an analytical interface condition (see chapter 6) 
which is highly sensitive to all heat fluxes except advection 
considered contributory to the thermal setting at the surface. 

A justification for the omission of advective heat flux has 
been offered elsewhere CPallmann et alii, 1973). 

In the subsequent paragraphs of the section, some of the 
characteristic output of the attuned F/0- and M/T-models will 
be presented as they pertimn to the dust-free case (February 16, 
1972). In choosing a M/T-generated temperature sounding, 
given in figure 6, as input to a F/0 -simulation instead of the 
IRIS-inverted T-sounding, we obtained again a high quality 
match with the measured IRIS data. The matching result is 
shown in figure 7 . For both curves , the pertinent surface 
temperature is T g = 261. 8K. 

A decomposition of the simulated spectral flux, welling 

up to the top of the atmosphere, has permitted us to illustrate 

the contributions of all atmospheric layers to that flux in 

the lower curve of figure 8. It is evident that the highest 

thermal level was achieved in the middle of both wings of the 

—1 

bulk CO^ absorption band centered at 667 cm * This atmos- 
pheric contribution to the spectral flux at the top was com- 
pared with the measured IRIS spectrum of February 16 in 
figure 9 . 

The simulated emission from the ground, subsequently 
attenuated by the atmosphere, has been related to the IRIS 




Fig. 7: Simulated flux and measured IRIS flux versus 

wave-number; the latter flux curve is marked 
by C*). F/0 input: T-sounding, date and 

location of Fig. 6. 
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Fig. 8: Simulated flux originating from ground and 

atmospheric layers Cupper curve) and strictly 
atmospheric contributions to flux at the top 
of Martian atmosphere (lower curve) versus 
wave-number; date and location: see Fig. 6. 
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Fig. 9: Atmospheric component flux (lower curve) as 

simulated, and measured IRIS flux (*) versus 
wave number; date and location: see Fig. 6. 


i 

f 
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Fig. 10: Simulated ground flux, subsequently attenuated 

by atmosphere (solid curve) and measured IRIS 
flux C*> versus wave-number; date and location: 
see Fig. 6. 
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data in figure 10. The area between the 2 curves is equal 
to that under the lower curve of the previous figure. Out- 
side the CO 2 -absorption band, the IR radiation is almost 
completely transmitted to the atmospheric top. This is a 
sensible result if one keeps in mind that the model does not 
accommodate photon-matter interactions with minor constituents 
and the spectral dependence of emissivity of silicate dust 
at the ground. 

In figure 11, the spectral flux transmissivity through 
the total depth of the atmosphere has been plotted. The de- 
cline of the curve below the value Ci.O) Is strictly de- 
pendent on atmospheric CO^. Particular attention Is called 
to the opacity from about 640 to 690 cra"^. 

In the remaining paragraphs of section 7.1, we will dis- 
cuss some typical outputs , generated by our time-dependent 
M/T model, either used further as an input into the F/0- 
simulation or simply presented in their own right. The former 
procedure permits us to elucidate more the dependence of the 
F/0 simulation on particular input data, whereas the latter 
may indicate the inherent capabilities of the M/T model. 

The program input parameters are given in Table 1. 

After a time -integration over more than 300 time-steps, 
the M/T model produced temperature soundings Ctautochrones ) 
for atmosphere and coil in balance for a 24-hour cycle 
(cyclical balance). The atmospheric sounding which evolved 
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Fig. 11: Spectral transmissivity of IR flux, through 

total depth of atmosphere, versus wave-number. 
Further specifications: Fig. 6, 
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TABLE 1 

FEBRUARY 16, 1972 

MARTIAN ATMOSPHERE RAD-COND-CONV HEAT TRANSFER SIMULATION 
PROGRAM PARAMETERS : 

INITIAL TIME 12.720 LOCAL MARTIAN SOLAR TIME 

NUMBER OF TIME STEPS 400 

TIMESTEP INTERVAL 10.000 MINUTES (MARTIAN) 

■4.000 MINUTES (MARTIAN) FROM 13.00 TO 17.00 
SURFACE ALBEDO .300 

SOIL DENSITY . 1.500 GR/CU-CM 

SOIL CONDUCTIVITY 2.300E 04 ERGS/ (SEC CM DEG) 

SOIL HEAT CAPACITY 1.005E ERGP/CCU-CM DEG) 

CO 2 CONDUCTIVITY 1.330E 03 ERGS/(SEC CM DEG) 

LATITUDE -41.40 DEG 

DECLINATION OF SUN -5.80 DEG 

DURATION OF SUNSHINE 760.8 MINUTES (MARTIAN) 

RATIO OF MEAN/ ATUAL. 

SOLAR DISTANCE 1.0002 

SURFACE PRESSURE 4.6 MB 

DUST AMPLITUDE 0 /CM 

ONE MARTIAN SOLAR DAY = 24 MARTIAN HOURS - 1440 MARTIAN MINUTES 
ONE MARTIAN MINUTE = 61.6491 SECONDS 
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Fig. 12. Cyclically balanced temperature sounding vs. 

height, in simulation for February 16, 1972, 
12.72 LMST, as input to the F/0 program, For 
additional parameters, see table 1. 
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for simulation day 3, 12.72 LMST as cyclically balanced has 
been plotted in figure 12. Three distinct characteristics 
appear in the temperature distribution. First, in the sur- 
face boundary layer of 12.5 m of thickness, the temperature 
drops by about 19 K with a balanced T » 256K. By com- 

5 

parison, the drop in the IRIS-inverted temperature sounding 
of the same local time is steeper* 

The second characteristic in figure 12 is the convec- 
tive layer which extends with its adiabatic lapse rate to 
a height of about 6 km. This sublayer is not discernible 
in the IRIS-inverted temperature curve. Thirdly, the upper- 
most branch of the simulated T-function in figure 12 has the 
typical shape of temperature distribution for a layer in 
which radiative transfer is the predominant heat transport 
mechanism. 

The dependence on wave-number of simulated spectral 
flux welling up at the top of the atmosphere, is illustrated 
in figure 13, The upper curve represents the contributions 
from ground and atmosphere duly attenuated by the intervening 
layers. The exclusive contribution to outward flux by the 
atmospheric strata is shown by the lower curve. The 
difference between the two curves at any wave-number indi- 
cates how much of the IR radiation, emitted by the ground 
surface, reaches the top of the atmosphere. 

Figure 14 mainly demonstrates the consequence of different 
surface temperatures implicit in the IRIS spectrum. (261. 810 
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hemisph/cm) vs. wave number in ClO/em) for 
February 16, 1972, 12„72 LMST* Upper curve: 
flux from ground and atmosphere combined; lower 
curve: exclusively atmospheric contribution to 
flux. For additional parameters, see table 1. 
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Fig, 14, F/Q simulated (solid curve) and IRIS-measured (*) 
up welling flux* in ergs/ ( see cnv 2 hemisph/cm) vs, 
wave number in (10/cm) for February 16, 1972, 
12.72 LHST. For additional parameters, see 
table 1* 
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^*ig. 15. F/0- simulated atmospheric contribution (lower 

curve) and IRIS-measured (*) upwelling flux 

in ergs/ (sec cm 2 hemisph/cm) vs. wave number 

in (10/cm) for February IB, 1972, 12.72 LMST 

for additional parameters, see table 1: 

T - 256K, 
s 
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and the cyclically balanced simulation spectrum (25610. In 
the center of the CO^ absorption band* the two curves co- 
incide. A comparison of the simulated atmospheric component 
flux with the measured IRIS flux spectrum is facilitated in 
figure 15. Between about 640 and 690 cm”' 1 ', the upwelling 
radiation is strictly originating from the atmosphere. No 
ground emission reaches the top of the atmosphere. 

As to the direct M/T model outputs for the dust-free 
condition of the southern-hemispheric midlatitudes, “he 
gradual evolution toward cyclical balance may be recognized 
In figures 16 and 17. In both figures, four-hourly atmos- 
pheric temperature soundings are presented. The sequence 
begins with the IRIS -inverted T-sounding at 12.72 LMST of 

February 16, 1975 and T - 259* 5K (IRR brightness temperature) 

s 

as input . 

Figure 16 shows as the last four-hourly sounding the 
curve at 12:00 LMST of day 2. Figure 17 is the continu- 
ation completing the second diurnal cycle at 12.72 LMST. It 
is evident that the curves of the second cycle stay closer 
together, particularly between 4 and 10 km. The surface 
temperature spread covers about 70K. The greatest thickness 
of the convective layer is reached at 12.72 LMST of day 3 
with 6 km. 

The first and last T-sounding may be compared in figure 
18. Since the sounding of day 3 is close to the 24-hour 
cyclical balance, we submit that it is more representative of 
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Fig, IB, M/T-simulated four-hourly atmospheric temperature soundings in degrees K vs. height 
in km for February 16, 1972, 12.72 LMST. Sequence of 24 hours begins with IRIS- 
in verted sounding and proceeds as indicated by markers Cupper right). For ad- 
ditional parameters, see table 1. 
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Fig. 17, Continuation of M/T-simulated four-hourly atmospheric temperature soundings in 

degrees K vs. height in km for February 18, 1972, 12,72 LWST. Sequence of 24 hours 
begins with IRIS-inverted sounding and proceeds as indicated by markers (upper 
right). For additional parameters, see table 1. 



Fig. 18. Initial and final soundings from the 48 hour M/T simulation, in degrees K vs. 

height in km, for February 13, 1972, 12.72 LMST . Markers (upper right) specify 
the times.,-. For additional parameters, see table 1. 
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the thermal stratification in the southern-hemispheric mid- 
latitudes near the end of the summer season* The largest 
differences of temperature for these two soundings are be- 
tween 4 and 25 km of height with up to 15K. As mentioned 
earlier, the temperature spread in the surface boundary layer, 
from 0 to 12.5 m, is about one third smaller for the M/T 
simulated sounding or about 20K. 

In figures 19 and 20, the four-hourly tautochrones (see 
Geiger, 1968) are presented for the period of 48 hours. The 
higher vertical resolution in scale reveals more detail. 

Steep near-surface inversions of up to 30K are discernible 
near sunrise. During the afternoon, the thickness of the con- 
vection layer clearly goes beyond 1,500 m. As to the sub- 
surface strata, below 60 cm the temperatures are assumed con- 
stant and equal to 222K. Figure 21 facilitates comparison be- 
tween the initial and the final tautochrone. One may be 
mindful of the limited vertical resolving power of the inte- 
gral inversion technique which was applied to measured IRIS 
spectra in order to obtain temperature soundings (Conrath, 
1974; see also chapter 5). 

7.2 ATTUNED MODEL OUTPUTS, MODERATE DUST-LOAD 

By inserting the height-dependent effective dust absorp- 
tion coefficient for December 24, 1971 (chapter 4) in the 
transmission function of both the F/0- and the M/T-simulation 
programs, we generated a variety of appropriate outputs which 
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Pig. 19. M/T-simulated four-hourly 

tautochrones in degrees K vs. 
height in decameters and depth 
in cm for February 16, 1972, 
12.72 LMST, Markers (upper 
right) indicate times. For 
additional parameters see 
table 1. 
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Initial and final tautochrones from 
the 48 hour M/T simulation, in degrees 
K vs. height in decameters and depth 
in cm, for February 16, 1972, 12:72 
LMST. Markers Cupper right) specify 
the times. For additional parameters, 
see table 1 , 
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Fig. 22. Cyclically balanced temperature sounding vs. 

height, in simulation for December 24, 1971, 
14.05 LMST, as input to the F/0 program. 

For additional parameters, see table 2. 
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were plotted via CALCOMP. The coefficient used to describe 
the effective absorption by dust was (1.2 x 10 - exp 

(-1.11 x 10 -9 .Z^)) cm”’ 1 '. For further detail, see section 4.3 
of this report. 

In the subsequent paragraphs, we will again present first 
some results generated by the F/ 0-program followed by out- 
puts pertinent to the M/T -model. The temperature sounding 
used as input for the F/ 0-simulation is given by figure 22. 

It is as before a cyclically balanced sounding at 14.05 LHST 

The surface temperature reads T = 264K; and, within the 

s 

surface boundary layer, a decline of 14K is discernible. 

There are two layers of convective overturn, one up to about 
6.5 km and the other from the top of the dust layer at 40 km 
down to about 27 km or twice as deep. At the top of the dust 
layer, the temperature is only 124K, a value gradually reached 
by the model In its trend toward a balance over a 24-hour 
cycle* 

The effective dust absorption coefficient is presented 
in figure 23. As stated earlier, the mathematical form of 
this function is Gaussian with a drop-off at 40 km (top of 
dust layer). In figure 24, the simulated flux is compared 
with the measured IRIS spectrum for this date. It is evident 
that the thermal level of the simulation curve is lower as 
a consequence of the reduced surface temperature in cyclical 
balance. Details of comparison within the CO^ absorption 
band will be offered in section 7.3 which covers the optically 
thicker dust condition. 
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Fig. 23. Effective dust absorption coefficient in 10” 
cm-i vs. height in km } for December 24, 1971. 
14.05 LMST. Top of the dust layer located at 
40 km. For additional parameters , see table 2 
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Fig. 24. F/0 -simulated (solid curve) and IRIS-measured (*) 

up welling flux, in ergs/ (sec cm 2 hemisph/cm) vs. 
wave number in (10/cm) for December 24, 1971, 
14.05 LMST. For additional parameters, see 
table 2. 
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hemisph/cm) vs. wave number in Cl 0/cm) for 
December 24, 19 71, 14.05 LMST, Upper curve: 
flux from ground and atmosphere combined; lower 
curve: exclusively atmospheric contribution to 
flux. For additional parameters, gee table 2. 
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In figure 25 , the simulated flux from ground and atmos- 
phere together with just the atmospheric contribution are 
shown. Outside the CO^ absorption band, centered at 667 cm" 1 , 
the involvement of the dust is clearly evident. A Planckian 
function of low temperature can be found to represent 
effectively the dust emissions. Again the highest thermal 
levels of the atmospheric emission upward may be recognized 
in the middle of the two wings of the CO^ band. The trans- 
mission capability, integrated over the entire depth of the 
atmosphere is shown in figure 26 as a function of wave number. 

In the remaining paragraphs of the section, typical out- 
puts of the M/T-simulation are to be presented. The program 
input parameters are given in Table 2 . The four hourly 
temperature soundings over a time period of 4-8 hours are 
illustrated in figures 27 and 28. Similar to the dust-free 
condition on February 16, 1972 (figures 16 and 17), a trend 
toward cyclical balance is discernible. The evolving temper- 
ature sounding was discussed above in conjunction with the 
F/ O-simulation of December 24, 1971. 

Four regimes may be distinguished in figure 28, i.e. , 
the planetary boundary layer with either nocturnal inversions or 
convection around local noon, the layer from 6.5 to 27 km 
with hydrostatic stability, the second convection layer be- 
tween 27 and 40 km, and finally a very stable layer topping 
the dust stratum. By comparing the thermal characteristics 
between 12 and 30 km with those for the dust-free situation 
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Fig. 26. Spectral transmissivity, for the entire depth of 
the atmosphere, vs. wave number in (10/cm) for 
December 24, 1971, 14.05 LMST. For additional 
parameters , see table 2 . 
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TABLE 2 

DECEMBER 24, 1971 

MARTIAN ATMOSPHERE RAD-COND-CONV HEAT TRANSFER SIMULATION 
PROGRAM PARAMETERS: 

INITIAL TIME 14*050 LOCAL MARTIAN SOLAR TIME 

NUMBER OF TIMESTEPS 300 

TIME STEP INTERVAL 10.000 MINUTES (MARTIAN) 

4.000 MINUTES (MARTIAN) FROM 13.00 to 17.00 
SURFACE ALBEDO .300 

SOIL DENSITY 1.500 GR/CU-CM 

SOIL CONDUCTIVITY 2.300E 04 ERGS/(SEC CM DEG) 

SOIL HEAT CAPACITY 1,0 05E ERGS/ (CU-CM ' DEG) 

CO 2 CONDUCTIVITY 1.330E 03 ERGS/(SEC CM DEG) 

LATITUDE -54,00 DEG 

DECLINATION OF SUN -16.20 DEG 
DURATION OF SUNSHINE 913.1 MINUTES (MARTIAN) 

RATIO OF MEAN/ATUAL. 

SOLAR DISTANCE 1.0483 

SURFACE PRESSURE 4.8 MB 

DUST AMPLITUDE 1.200F-07 /CM 

ONE MARTIAN SOLAR DAY =24 MARTIAN HOURS - 1440 MARTIAN MINUTES 
ONE MARTIAN MINUTE = 61.6586 SECONDS 
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Fig. 2 7. M/T-simulated four-hourly atmospheric temperature soundings in degrees K vs. height 
in km for December 24, 1971, 14.05 LMST. Sequence of 24 hours begins with IRIS- 
inverted sounding and proceeds as indicated by markers (upper right). For ad- 
ditional parameters", see table 2. 
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Fig. 28, Con t in u at ion of M/T -simulated four-hourly atmospheric temperature soundings in 

degrees K vs. height in km for December 24 3 1971, 14,05 LMST, Sequence of 24 hours 
begins with IRIS-in verted sounding and proceeds as indicated by markers Cupper 
right). For additional parameters, see table 2. 
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to 2 OK due to the presence of dust. The intense cooling near 
the dust top which is related to the discontinuity of the 
radiatively active material concentrations, may be partially 
caused by the particular modeling procedure about that level. 

In figures 29 and 30 , the four-hourly tautachrones for 
December 24 are plotted. Because of the higher vertical 
resolution, detailed structure characteristics of the nocturnal 
temperature inversions or convective stratification appear. 

The thermal conduction in the soil depends on the values of 
various related parameters, given in TABLE 2. Initial and 
final soundings of this time-dependent run are fairly similar 
because a simulated sounding was used as input, 

7.3 ATTUNED MODEL OUTPUTS , OPTICALLY-THICK DUST-LOAD 

The day of November 22, 1971 was chosen to represent the 
condition of the opticaliy-thick dust-load in the southern- 
hemispheric, mid-latitudes , The particular value chosen for 
the height-dependent coefficient of effective absorption by 
dust particulates was (3.0 x 10“^. exp (-1.11 x 10“ 9 .Z 2 )) cm” 3 *. 

It determined the effective dust transmission, as discussed 
in chpt . 4 of this report. Its height-dependence is illustra- 
ted in figure 31. A sharp drop-off is discernible at 40 km 
marking : the adopted top of the dust-layer. 

The input sounding for the F/0 -simulation is shown in 
figure 32, This temperature sounding is the result of more . . 
than 300 time-integrations with the M/T-progranu Cyclical 
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Fig, 29. M/T- simulated four-hourly 

tautochrones in degrees K vs. 
height in decameters and depth 
in cm for December 24 , 1971, 
14.05 LMST, Markers Cupper 
right) indicate times. For 
additional parameters see 
table 2 . 
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Fig. 31, Cyclically balanced temperature sounding vs. 

height, in simulation for November 22, 1971, 
9.52 LMST, as input -to the F/0 program. For 
additional parameters, see table 3. 
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Fig, 32 * Effective dust absorption coefficient in 10 

cm"i vs* -height in km, for November 22, 1971, 
9.52 LMST, Top of the dust layer located at 
40 km* Fox* additional parameters, see table 3. 
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No vember 22 , 19 71, 9.52-LMST. Upper* curve: 
flux from ground and atmosphere combined; lower 
curve: exclusively atmospheric contribution to 
flux. For additional parameters, see table 3. 


- 102 - 

balance over 24 hours is satisfactorily approximated. More 
details of its structure will be presented in conjunction with 
the discussion on the 4-hour ly soundings over a 2-day period, 
near the end of this section. In figure 33, the F/0- simulated 
flux, emitted by ground /atmosphere and subsequently transmit- 
ted to its top Cupper curve), as well as the strictly atmos- 
pheric contribution (lower curve) are depicted as a function 
of wave-number. The percentual fraction of the atmospheric 
spectral flux appeared to be doubled when compared with that 
of December 24, 1971. The difference between the 2 curves, 
shown in figure 33, is a measure of the ground emission which 
arrives at the top of the atmosphere. Clearly, near 667 cm”^, 
no surface emanation reaches the top. The dust-laden C0 2 ~ 
atmosphere has effectively absorbed all photons near this wave 
number. 

A comparison between the simulated and measured upwelling 
fluxes is offered in figure 34. The difference in the levels 
of the curves outside the C02~bands is chiefly related to the 
difference in surface temperatures T g ." The IRIS-T s = 238K is 
more than 20K lower than that of the cyclically balanced soun- 
ding. As will be seen in figure 38, the IRIS -inverted tempe- 
rature sounding is not in cyclical balance. Within the C0 2 ~ 
bands around 667 cm“^, the simulated flux spectrum shows too 
low a temperature level. Various interpretations for this 
phenomenon can be offered. One relates it to the omission 
of the Doppler shape factor (see sect. 3*3). This shape fac- 
tor makes a better allowance for gappiness between C0 2 absor- 
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Fig. 34 


F/0 simulated (solid curve) and IRIS -measured (*) 
upwelling flux, in ergs/ (sec cm^ hemisph/cm) vs. 
wave number in (10/cm) for November 22, 1971, 

9.52 LHST. For additional parameters, see 
table 3. 


-104- 



Fig. 35. Spectral transmissivity;, for the entire depth of 
the atmosphere, vs. wave number in (10/cm) for 
November 22, 1971, 9.52 LMST. For additional 
parameters, see table 3. 
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ption lines. Wider spectral Raps permit a larger number 
of IR-photons from lower, therefore warmer, regions to leave 
the atmosphere. This intensified transmission would lead to a 
lift of the simulated flux curve within the C0 9 bands. 

Another interpretation makes reference to a study of ter- 
restrial mesopheric dust granule temperatures by Fiocco et 
alii (1975) who deduced dust temperatures (up to 100K) higher 
than the ambient temperature. If a similar condition exists 
on Mars we would be dealing with hotter particulate emitters 
embedded in a colder environment. This would generate an 
enhanced emission measurable as a higher upward flux value. 

The spectral transmissivity for the total depth of the 
atmosphere (50 km) is depicted in figure 35 as a function of 
wave-number. A total opacity from. 640 to 690 cm"'^ is again 
discernible. Outside the CO 2 *- bands, about 42 % of the I li- 
ra di at ion from the ground arrives at the top of the atmos- 

„ n 

phere. Minor vacillations near 1,000 and 1,920 cm stem 

from CO 2 -absorption. 

In the remaining paragraphs of section 7.3, M/T-simulated 
time-dependent outputs will be discussed. The adopted values 
of the program parameters are shown in Table 3. The time evo- 
lution of the atmospheric temperature stratification over 48 
hours has been plotted in figures 36 and 37. The sequence of 
4-hour ly temperature soundings starts with the IRIS-in- 
verted T-sounding whose surface temperature was set to = 

238K. Its temperatures between 15 and 30 km remain the lowest 
over the entire period of time integration. This suggests 
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TABLE 3 

NOVEMBER 22, 1971 

MARTIAN ATMOSPHERE RAD-COND-CONV HEAT TRANSFER SIMULATION 
PROGRAM PARAMETERS : 

INITIAL TIME 9.520 LOCAL MARTIAN SOLAR TIME 

NUMBER OF TIMESTEPS 300 

TIMESTEP INTERVAL 10.000 MINUTES (MARTIAN) 

if. 000 MINUTES (MARTIAN) FROM 13.50 TO 16.50 
SURFACE ALBEDO .300 

SOIL DENSITY 1.500 GR/CU-CM 

SOIL CONDUCTIVITY 2.300E 04 ERGS/ (SEC CM DEG) 

SOIL HEAT CAPACITY 1. 005E 07 ERGS/(CU-CM DEG) 

C02 CONDUCTIVITY 1.330E 43 ERGS/ (SEC CM DEG) 

LATITUDE -40.10 DEG : 

DECLINATION OF SUN -21.90 DEG 

DURATION OF SUNSHINE 878 .1 MINUTES (MARTIAN) 

RATIO OF MEAN/ACTUAL. 

SOLAR DISTANCE 1.0737 

SURFACE PRESSURE 5.3 MB 

DUST AMPLITUDE 3.000E-07 /CM 

ONE MARTIAN SOLAR DAY = 24 MARTIAN HOURS = 1440 MARTIAN MINUTES 
ONE MARTIAN MINUTE = 61.6625 SECONDS 
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Fig. 36. M/T-simulated four-hourly atmospheric temperature soundings in degree K vs, height 
in km for November 22, 1971, 9.52 LMST . Sequence of 24 hours begins with IRIS- 
in verted sounding and proceeds as indicated by markers (upper right). For ad- 
ditional parameters, see table 3. 
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Fig. 37, Continuation of M/T-simulated four-hourly atmospheric temperature soundings in 

degrees K vs. height in km for November 22, 1971, 9.52 LMST, Sequence of 24 hours 
begins with IRIS-inverted sounding and proceeds as indicated by markers (upper 
right). For additional parameters, see table 3. 
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that radiative heat transfer mechanisms, maintain in this 
layer temperature levels incresed by up to 30K. Above and 
below this layer,, additional heat transport accurs by 
means of convective overturn, characterized by the adiabatic 
lapse rate of 5.2K/km (straight lines slanting from the upper 
left to the lower right) „ 

Again, temperature inversions are discernible in the lo- 
west atmosphere around sunrise, and above the dust-layer top 
(at HO km) with variable strength over the diurnal cycle. 

In figure 37, the approach toward the cyclical balance may 
be recognized. Diurnal temperature ranges are the following: 

- at the surf ace , about 7 IK ; 

- in the lowest 2 km, about 4 OK; 

- at the top of the planeatry boundary layer, about 15K; 

- near 30 km, about 22K; and 

- in the uppermost atmospheric layer above the dust, 
about 6K, 

If a comparison is made between these simulated tempera- 
ture characteristics and the S-band occultation soundings of 
temperature , it is apparent that similar basic features exist 
in the occultation data, for instance, the S.H. midlatitudi- 
nal occultation temperature profile of revolution 17 (Kliore 
et alii, 1972; fig. VI- 3) shows a temperature minimum of about 
160K at Martian radius 3,418 km, with a quasi-adiabatic stra- 
tification, several km deep, below this level. Furthermore, 
the lapse rate of the occultation profile in the lowest kilo- 
meters appears to be quasi-adiabatic, whereas the rate of 
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Fig. 38. Initial and final soundings from the 48 hour K/T simulation, in degrees K vs . 

height in km, for November 22, 1971, 9.52 LMST. Markers Cupper right) specify 
the times. For additional parameters, see table 3. 
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temperature change in the intermediate layers is just a small 
fraction of the adiabatic lapse rate. 

This coincidence of some features in the simulation and 
occultation soundings appears to be repeated for profiles of 
revolution 1, 3 , 7 , and 15 (Kliore’s fig. VI-2). Although it 
is recognized that the occultation data are noisier, in prin- 
cipal they have a higher vertical resolution than IRIS-inver- 
ted temperature soundings (see chapter 5 of this report):* 

In figure 38, the approach toward the cyclical balance 
is illustrated by showing the initial (IRIS-in verted) and the 
final (simulation) temperature soundings for the period of 48 
hours. It is obvious that a rather large discrepancy exists 
between 10 and 35 km. This is a consequence of the simulated 
heating by dust, which is incorporated into the M/T-model. 

If we assumed that the minimum in the upper portion of the 
occultation temperature profile is at about 40 km, the pro- 
file would follow more closely the final sounding. Conversely, 
if the height of the dust layer had been simulated at 32 km, 
the M/T-generated final sounding would have appeared similar 
to the smoothed occultation sounding. 

In figures 39 and 40, the M/T simulated four-hourly 
tautochrones are depicted for November 22_, 19 71. The initial 
sounding at 9.52 LMST is the IRIS-in verted temperature dis- 
tribution. The higher spatial resolution of the height-scale 
permits us to discern more easily the near-sunrise thermal in- 
versions up to 1,500 m and the near-noon adiabatic 
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Fig. 40, Continuation of M/T-simulated 
four-hourly tautochrones in 
degrees K vs. height in 
decameters and depth in cm for 
November 22, 1971, 9.52 LMST. 
Markers Cupper right) indicate 
times. For additional param- 
eters see table 3. 
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stratification due to convection. Again, the second simu- 
lation day (Fig. 403 illustrates the approach toward cyclical 
balance with the last sounding at 9.52 LMGT of simulation- 
day 3 being taken to be sufficiently balanced. 

The sub-surface temperature variations of simulation 
day 2 and 3 clearly indicate that changes comparable in 
magnitude to the range of T s are limited to the first 10 to 
20 cm below the surface. This feature depends on the numeri- 
cal values, chosen for the thermal parameters of the silicate- 
soil (Table 3). For very steep subsurface temperature grad- 
ients, the flux of thermal conduction assumes just below the 
surface rather large values which are a major fraction of 
those for the noon insolation. This indicates the heat 
storage and discharge capability of the Martian soil over the 
24-hour cycle. The area embraced by the two tautochrones 
with highest and lowest temperatures may be taken as a measure 
of the oscillatory heat transports of the soil. 

For easier comparison, the initial and final tautochrones 
are shown in figure 41. The trend toward cyclical balance in- 
creased both the atmospheric and soil temperatures bv more 
than 10K almost everywhere. The tautochrones of figure 41, 
taken together with the soundings of figure 38y delineate the 
-temperature distributions of the ground-atmosphere for 
cyclically balanced and unbalanced conditions. 
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Fig. 41. Initial and final t auto ah rones from 

the 48 hour M/T simulation, in degrees 
K vs. height In decameters and depth 
in cm, for November 22, 1971, 9.52 LMST 
Markers Cupper right) specify the 
times. For additional parameters, see 
table 3. 
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8.0 SUMMARY OF FINDINGS 

The attunement of the refined model to IRIS, IRR 
and S-band o.ccultation data facilitated the creation of 
a broad spectrum of outputs (chapter 7) which has been 
serving as an information pool for interpretation. With 
respect to the atmospheric dust-load, it permitted us 
to generate inferences into the representativeness of 
' IRIS -inverted temperature soundings, the significance of 
IR-radiometric brightness temperatures, and the useful- 
ness of S-band occultation data. 

The underlying principle for our time-dependent heat 
transport simulation was that atmospheric and soil tempera- 
ture fields are the result of free exchange of heat 
through radiative, convective and conductive transfers. 

In this context, advection was not explicitly modeled into 
our program. If given sufficient time, the temperature 
fields generated by the R.C.C.T, model will tend to approxi 
mate a balance over the 24-hour cycle. 

In the following paragraphs, some of the more basic 
findings of the data-attuned simulation effort are given: 

1) In the dust-free atmosphere, a very close attune.- 
ment of the combined M/T-F/ 0-model is attainable. 
It is inferred that the cyclically balanced 
temperature sounding, generated for the dust-free 
condition, is more representative of the actual 
thermal field than the IRIS-inverted sounding if 
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there is no advective interference with the cyclical 
balance. Conversely, one may conclude, in the 
presence of advective interferences, deviations from 
the cyclically balanced sounding are likely to 
occur. 

2) The cyclically balanced temperature sounding basically 
reveals, 3 regimes in the vertical structure of the 
Martian atmosphere, i.e., (a) a surface boundary layer 
in which diurnal temperature variations of up to 
about 7 OK (at the ground) occur; (b) a planetary 
boundary layer whose top is determined at about 5 km 
by intense radiative-convective heat transfer; and 

(c) the "free 11 atmosphere with predominantly radiative 
heat transfer in which diurnal temperature variations 
are less than 10K. 

3) For moderate dust-loads (such as on December 24, 1971) 
the attunement of the combined M/T-F/ 0-model is fairly 
close. It is inferred that the cyclically balanced 
temperature sounding, generated for this case appears 
in its detail meteorologically more realistic than 
the IRIS -inverted sounding if advective interference 
with the cyclical balance is discarded. Conversely, 
we may conclude that advection could cause deviation 
from the cyclical balance. 

4) For a dust layer of moderate concentration, the 
cyclically balanced temperature sounding shows 4 basic 
thermal regimes in the vertical structure of the 
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Martian atmosphere, i.e., (a) a surface boundary layer 
with a diurnal temperature variation of up to 68K (at 
the ground) ; (b) two convective-radiative transfer 
layers, one from 0 to 6 km and the other from the top 
of the dust layer (40 km) down to 27 km; (c) a layer 
from 6 to 27 km in which radiative transfer is the con- 
trolling transport mechanism; and (d) a rather intense 
(40K) inversion layer topping the dust layer. 

5) For optically- thick dust-loads (November, 1971), the 
attunement of the combined M/T-F/O-model is fair 
enough to serve as a basis from which to infer. We 

find again that because of its higher vertical resolution 
the cyclically balanced temperature sounding is meteoro- 
logically more representative than the IRIS-inverted 
sounding unless advective heat transport interfere. Con- 
versely, it may be concluded that, in the presence of 
effective advection, deviations from the balanced sounding 
are to be expected. 

6) For an optically-thick dust layer with a diumally 
balanced thermal cycle, the generated sounding indicates 
again 4 basic regimes in the vertical atmospheric struc- 
ture, i.e., (a) a surface boundary layer with a diurnal 
temperature variation of up to 72K (at the ground); (b) 
two convective-radiative layers, one from 0-8 km and the 
other from the top of the dust layer (40 km) down to 
about 32 km; ( c) a layer from 8- 32 km with radiative 
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transfer dominating the heat exchange and causing diurnal 
temperature variations of at least 15K; (d) a rather 
steep inversion layer of between 6 and 18 K temperature 
increase topping the dust layer. 

7) The Martian dust particulates appear to be rather effec- 
tive absorbers with a likely complex index of refraction 
m = 1.55 - i (0.0465 ± 0,0155) for the visible range. 

8) The diurnal history of the cyclically balanced tempera- 
ture sounding does permit us to re-interpret the IR- radio- 
metric (IRR) brightness temperature. As an example j on 
November 22, 1971, 9.52 LMST , the averaged IRR-temperature 
was 234. 7K whereas the balanced surface temperature showed 
T = 259K. In our sounding, the value 234. 7K was reached 
at an elevation of about 800 m. 

9) In spite of its inherent noisiness, the S-band occultation 
data confirm some of the features in the vertical tempera- 
ture structure, as generated by our simulation. Conversely, 
our cyclically balanced T-soundings may be used to adjust 
occultation temperature soundings. 
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